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In the ILD detector, HCAL is inside the coil. The HCAL is about 5 interaction length 
thick, and ECAL is about 1 interaction length thick. For example, an eighty GeV 
hadron, it is only 95% energy will be contained in ECAL plus HCAL. We need a 
topological reconstruction of the leakage. This has been studied with experimental data 
collected using the CALICE prototypes, during the 2007 test beam campaign at CERN. 
The complete setup of the experiment consisted of a silicon-tungsten electromagnetic 
calorimeter, an analog scintillator-steel hadron calorimeter and a scintillator-steel tail 
catcher. Events collected using pion beams in the energy range 8-100 GeV are selected 
and compared to the Monte Carlo simulations. While the leakage from the full setup 
is negligible, when removing the tail catcher information either partly or completely 
the energy loss becomes significant and affects the performance. The average measured 
energy decreases below the expected beam energy and the resolution deteriorates. A 
correction to the leakage was implemented for pions. The results obtained show that the 
correction is powerful in restoring the mean value of the measured energy distributions 
back to the expected beam energy, with an accuracy at the 1-2% level over the whole 
energy range. The relative improvement on the resolution is about 25% at 80 GeV, 
decreasing at lower energies together with the impact of the leakage. 

1 The AHCAL prototype detector 

The goal of the CALICE experimental pro- 
gram is to establish novel technologies for 
calorimetry in future collider experiments 
and to record electromagnetic and hadronic 
shower data with unprecedented three di- 
mensional spatial resolution for the valida- 
tion of simulation codes and for the test 
and development of reconstruction algo- 
rithms. Such highly granular calorimeters 
are necessary to achieve an unprecedented 
jet energy resolution at the International 
Linear Collider |1 using particle flow al- 
gorithms. [2] The CALICE analog hadron 
calorimeter (AHCAL) [3] prototype is a 38 
layers sampling calorimeter, which is built 
out of scintillator tiles with sizes ranging 
from 30 x 30 mm 2 in the core of the de- Figure 1; Picture of one AHCAL active layer 
tector to 120 x 120 mm 2 . The light in each 
scintillator cell is collected by a wavelength 

shifting fiber, which is coupled to a silicon photomultiplier (SiPM) SiPM. The SiPMs are pro- 
duced by the MEPhl/PULSAR group. [4] They have a photo-sensitive area of 1.1 x 1.1 mm 2 
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containing 1156 pixels with a size of 32 x 32 /im 2 . In total, the calorimeter has 7608 channels. 
A built-in LED calibration system with UV LEDs, is coupled to each cell by clear fibers, 
and equipped with PIN diodes to monitor the LED light intensity. The performance of the 
AHCAL was validated with positrons at various energies [3] . Good linearity and satisfactory 
agreement with simulation up to an energy of 50 GeV has been observed. [5j|6] 

2 Identification of track segments 

One fine granular AHCAL prototype layer is shown in figure [T] The high granularity of 
the active layers in the hadronic calorimeter and the cell-by-cell readout gives the CALICE 
detectors unprecedented imaging capabilities. 

Thanks to the high granularity, the 
development of hadronic showers can be 
studied in great detail. The response of 
the AHCAL to hadronic showers are in- 
vestigated. This is exploited to study 
the topology of hadronic events in detail. 
Track segments created by charged particles 
produced within the hadronic shower can 
be identified, provided the particles travel 
an appreciable distance before interacting 
again and are separated from other activity 
in the detector. 

Figure [2] demonstrated that three tracks 
were found within one 25 GeV 7r event. The 
properties of these tracks are sensitive to the 
substructure of the hadronic shower, and 
can thus serve as a powerful probe for hadronic shower models. The track segments identified 
in hadronic showers have a high quality, and are also suitable for detector calibrations via 
the extraction of the most probable value of the energy loss in each cell along the track. [7] 

3 Longitudinal shower profile 

Hadronic shower development is based largely on nuclear interactions, and therefore hadrons 
can travel a significant distance before a shower develops. The average distance a hadron 
travels before a nuclear interaction occurs, which typically starts the cascade of the shower, 
is called nuclear interaction length Xi n t- 

The high granularity of the calorimeter provide the possibility to determine the first hard 
interaction of the showers, it is possible to reconstruct the development of the showers relative 
to their starting point, instead of considering the development relative to the calorimeter 
front face. When the shower starting point for each event is known, it is possible to shift the 
single event profile to the same point. The plot in Figure [3] shows the longitudinal profile for 
80 GeV pion showers. These two ways to reconstruct the longitudinal profile are compared. 
The open symbols indicate that the longitudinal shower profile becomes much shorter after 
correcting event-by-event for the variation of the shower starting point. The layer-by-layer 
effects are clearly washed out and the profile development appears to be smooth. 




Figure 2: Identification of minimum-ionizing 
track segments within hadronic showers. 
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Figure 3: longitudinal profile in the AHCAL only, drawn with respect to the first hard 
interaction of the showers or w.r.t. the calorimeter front face. Showers starting in the first 
layer of the AHCAL have been rejected, therefore the longitudinal profile relative to the 
shower starting point does not reach the very last bin. 



The measurement of the position of the first hadron interaction allows to disentangle 
the fluctuations of the position where the shower statrts from those intrinsic in the shower 
process. 



4 Leakage correction 
4.1 The leakage 

In figure |4j the distribution on the left plot shows the total energy measured using the 
full calorimeter, while the distribution on the right plot shows the energy measured when 
excluding the TCMT information. To quantify the energy resolution, the measure RMS90 
has been used, defined in [2] as the root mean square of the smallest interval containing 90% 
of the the distribution. The degradation is clearly visible on the distribution when excluding 
the TCMT information. It shows both in terms of mean and RMS90 in the distribution on 
the right plot. A topological reconstruction of the leakage has been studied and implemented 
to correct the leakage. 



4.2 Observables sensitive to the leakage 

Since the information of the TCMT is excluded, measured energy will indicate the energy 
measured by SiW-ECAL+ AHCAL, where not otherwise specified. Two observables have 
been studied which are sensitive to the leakage. They are shown in figure [5j The left 
distribution shows the correlation between the shower starting point and the ratio between 
the measured energy and the beam energy. When the shower starts at the beginning of the 
AHCAL, the measured energy is close to the expected beam energy for most of the events. 
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Figure 4: Left: total energy measured using the full calorimeter. Right: energy measured 
by SiW-ECAL+AHCAL only. 



The leaking energy from the AHCAL, is usually negligible. On the contrary, for showers 
starting late in the AHCAL, the fraction of the beam energy which is measured becomes 
smaller. That means only partial energy deposited in the AHCAL. And the effect of the 
leakage is non-negligible anymore. As already mentioned, due to event-to-event fluctuations 
typical for hadronic showers, the correlation between the leakage and the shower starting 
point is spoiled, resulting in a broad distribution. 

In order to recover the incident particle energy from the leakage, the distribution on the 
left in figure [5] shows the average correction to be applied to the measured energy, depending 
on the layer where the first hard interaction occurred. 

The right distribution in figure [5] shows the end-fraction versus the ratio between the 
measured energy and the beam energy, for an 80 GeV run. The end-fraction is defined as the 
fraction of energy deposited in the last four layers of the AHCAL divided by the total shower 
energy measured by SiW-ECAL+AHCAL. A small fraction of energy deposited at the end 
of the AHCAL indicates that the shower is in general almost concluded before reaching the 
end of the calorimeter. In fact, the ratio between the measured energy and the beam energy 
is mainly close to one, for those events that have a small end-fraction. On the contrary, a 
large end- fraction indicates a shower still in development, that could leak into the TCMT 
with a high probability. As a consequence, the ratio between the measured and the beam 
energy decreases when the end-fraction increases, due to increasing leakage. 



4.3 A correction for the Leakage 

A correction method has been investigated and implemented within the CALICE collabo- 
ration by Ivan Marchesini for the leakage. M The correlation between the leakage and the 
shower starting point, the end-fraction has been used together in order to develop a real- 
istic correction for the leakage. The shower starting point is calculated as the layer. The 
end-fraction is a continuous variable and can be arbitrarily binned. The chosen binning 
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Figure 5: Left: shower starting point versus the ratio between the measured energy in SiW- 
ECAL+AHCAL and the beam energy. Right: end-fraction versus the ratio between the 
energy measured in SiW- ECAL+AHCAL and the beam energy. The cluster of events for 
an end-fraction of approximately one is due to the fact that events with a shower starting 
point in layers 35 or 36 have an end-fraction exactly equal to one according to the end- 
fraction definition. 

is [0, 0.01, 0.025, 0.05, 0.1, 0.15, 0.2, 0.25, 0.3, 0.35, 0.4, 0.45, 0.5, 0.55, 0.6, 0.65, 0.7, 
0.75, 0.8, 0.85, 0.9, 0.95, 1]. The sum of measured energy in SiW-ECAL and AHCAL has 
been used to replace the beam energy information. In this case realistic means that the 
correction should not be tuned using the beam energy information. The energy can be cor- 
rected event-by-event depending on the measured energy, the shower staring point and the 
end-fraction. 

4.4 Correction performance 

The correction was first implemented in a Monte Carlo study. Two possible applications 
to the data have then been considered, obtaining the correction to the data either from a 
Monte Carlo template or from independent data runs. 

The performance of the correction is summarized on the top row in figure [6] for the 
Mont Carlo studies. The correction is capable of identifying the energy of the run with a 
precision better than 0.5%. The RMS90/Mean improves by up to 30% at high energies, 
while decreasing towards lower energies together with the impact of the leakage. 

The performance of the correction on the test beam data is shown on the bottom row 
in figure [6] Pion data collected using the CALICE prototypes during the 2007 test beam 
campaign at CERN were analyzed. When applying the correction to data, the looking up 
table can be obtained either using an independent set of data or Monte Carlo simulations. 
In the first option data from runs with energies of 8 GeV, 12 GeV, 15 GeV, 25 GeV, 35 GeV, 
45 GeV, 80 GeV and 100 GeV have been used. The correction has then been applied to 
events from independent runs with energies of 10 GeV, 20 GeV, 25 GeV, 35 GeV, 45 GeV, 
60 GeV and 80 GeV. The events from the three runs with energies of 10 GeV, 20 GeV and 
60 GeV allow the performance of the correction to be tested also at energies not used to 
build the looking up table. 

The results are shown in figure [6] The mean energy measured by SiW-ECAL+AHCAL 
is shifted to the expected beam energy, correcting for the average leakage, with an accuracy 
better than 1%, which is comparable to the precision given by the additional information 
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Figure 6: Top: Final results of the leakage correction for Monte Carlo. Left: the relative 
deviation of the mean value of the corrected energy distributions from the beam energy 
versus the beam energy. Center: mean value of the measured energy distributions before 
and after the correction. Right: relative improvement of RMS90/Mean versus the beam 
energy. Bottom: Final results of the leakage correction applied to data runs. Left: the 
relative deviation of the mean value of the corrected energy distributions from the beam 
energy versus the beam energy. Center: mean value of the measured energy distributions 
before and after the correction. Right: relative improvement of RMS90/Mean versus the 
beam energy. The templates for the correction have been obtained from independent data 
runs. 
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of the TCMT. The RMS90/Mean improvement is up to 25% at 80 GeV while decreasing 
towards lower energies together with the impact of the leakage. 

4.5 Conclusion 

The correction was applied to selected pion data and it was tuned cither using independent 
data runs or Monte Carlo simulations. The obtained results are better in case the correction 
is built from data, due to discrepancies between data and simulations. The correction 
reduces significantly the impact of the leakage. The mean of the energy distributions after 
the correction is centered around the expected beam energy with a 1% or 2% accuracy, 
depending on whether the corrections are built from data or from Monte Carlo simulations. 
Such a precision is approximately the same given by the additional use of the TCMT. The 
resolution is improved by up to 25% at high energies, where the impact of the leakage is 
more pronounced. The improvement decreases towards lower energies together with the 
impact of the leakage. 
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